fn,S>. &L./* . ,L 

E 


CONTRACTOR-PREPARED FINAL REPORT 
CONTRACT NASB-36 153 


Group-Kinetic Theory and Modeling of Atmospheric Turbulence 


Period: December 1, 19BA - November 30, 1987 


Prepared for George C. Marshall Space Flight Center 
Marsha 1 1 Space Flight Center , Alabama 35818 


- "7L -* r- THEORY A nD 

,, t * . -k.rr rinal 

' S ' 1G Ir«nSI>H1=RlC TlIRBUttNC-^ (0S „ A) 

NG OF ATrt ,o A - jO Nov. 1 C^Ct ORft 

l one. 1 


Unci as 

0 7 2 2 T 2 6 


CONTRACTOR: Universities Space Research Association 
Headquarters Office, American City Building Suite SIS, 
Columbia, Maryland SIO^A 


Report prepared by: Universities Space Research Association 
Boulder Program Office, P.O. Box 391, Boulder, Colorado 80306 


FINAL REPORT FOR CONTRACT NASB-36153 
"GROUP-KINETIC THEORY AND MODELING OF ATMOSPHERIC TURBULENCE 


Attached is the Final Report prepared by the Research Specialist* Dr . 
C. M. Tchen, who performed the research sponsored under Contract 
NAS8-36153. In this Report, Dr. Tchen documents the results of his 
three-year effort, as called for in Attachment J-S of the Contract. 

Also attached are two additional documents which describe Tchen’s 
analysis method in detail and apply it to Rossby Nave Turbulence. 


Respectfully submitted, 


M. H. Davis, Ph . D . 

Program Director for USRA 
August 8, 1989 


FINAL REPORT 


GROUP-KINETIC THEORY AND MODELING OF ATMOSPHERIC 

TURBULENCE 


By C. M. TCHEN 


NASA Contract No. NAS 8-36153 
Three years from November 30, 1984. 
June 1988 


Monitored by: Dr. G. Fichtl 

George C. Marshall Space Flight Center 
Marshall Space Flight Center, AL 35812 

Contractor: Universities Space Research Association 

American City Building, Suite 311 
Columbia, MD 21044 
. Dr. M. H. Davis 


Principal Investigator: 

Dr. C. M. Tchen 
900 West 190th Street 
New York, NY 10040 
(212) 781-4111 


4 


i 


GROUP-KINETIC THEORY AND MODELING OF TURBULENCE 


C. M. TCHEN 

1. Research objectives and method of treatment 

A group-kinetic method is developed for analysing eddy transport 
properties and relaxation to equilibrium. The purpose is to derive 
the spectral structure of turbulence in incompressible and compres- 
sible media. Of particular interest are: direct and inverse 

cascade, boundary layer turbulence, Rossby wave turbulence, two- 
phase turbulence; compressible turbulence, and soliton turbulence. 
Soliton turbulence can be found in large-scale turbulence, tur- 
bulence connected with surface gravity waves and nonlinear propaga- 
tion of acoustical and optical waves. 

Since the group-kinetic method is basic to our statistical 
theory of turbulence, two new Technical Reports are enclosed here: 

#1. Group-kinetic method of turbulence 

#2. Group-kinetic theory of Rossby wave turbulence. 

2. Group-kinetic theory of turbulence 

The most appealing method in the last three decades for treating 
turbulence statistically has been the ''renormalization perturbation 
expansion" and the "kinetic method of turbulence". The former has 
enabled Kraichnan (1959, 1977) to develop his DIA method and its 
generalization. A systematic method of developing the expansion in 
a general form was given by Martin, Siggia and Rose (1972). By the 
difficulty of convergent summation and by the indefinite complings 
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among modes without explicit representation of physical processes, 
the perturbation method has not succeeded in analytically deriving 
the eddy coefficients (eddy viscosity and eddy amplification) that 
can discriminate between the direct and inverse cascades. These 
difficulties become more evident when turbulence is characterized 
by more than one parameter such as shear turbulence, convective tur- 
bulence and Rossby wave turbulence. 

In fact, the molecular viscosity as a transport property is 
known to be of kinetic origin such that it has to be calculated 
from the Boltzmann equation (Chapman and Cowling, 1970) . If by 
analogy, the eddy viscosity has to be derived from the "collision 
operator", the kinetic equation of turbulence has first to be for- 
mulated. Any attempt that leads to many-point distributions 
(Bogoliubov, 1962; Struminskii, 1985) has not shown progress. 

The determination of the spectral solution of the Navier-Stokes 
equation by an analytic statistical method is important not only 
for reproducing the Kolmogoroff spectrum and evaluating the Kol- 
mogoroff constant, but also for understanding the explicit tran- 
sport proce'sses and transport coefficients. This examination will 
help in assessing the applicability of the statistical method to 
other problems of turbulence. 

For the description of the microdynamic state of turbulence, we 
let the pressure fluctuations govern the elementary interaction be- 
tween fluid elements, by writing the Navier-Stokes equation in the 
form 

(3^ + u’V -vV 2 )u = E, V*u = 0 (2.1) 
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The kinetic equivalent is called the caster equation 


(3 t + L)f (t,x,v) = 0, 


( 2 . 2 ) 


by f(t,x,v) = 6Cv-u(t,x) D . Here u is fluid velocity, 

E = - ~^P is the gradient of pressure p per unit mass-density p, 
and v is molecular viscosity. The velocity distribution is f in 
the phase-space t,x,v . The differential operator is 


vV-vV 2 +E*3, 3 = 3/3v 


An instantaneously fluctuating func- 


- ' — v 

tion is denoted by ()=() + (). and can be decomposed into an 

- _ *\, 

ensemble average ( ) = < > and a fluctuation ( ) . The scaling 

operator A and ^ can be used. 

We develop a group-kinetic method that incorporates the tran- 
sport processes of spectral evolution, eddy viscosity and relaxa- 
tion as macrogroup, microgroup and subgroup into the multi-scale 
distribution functions f°, f', f = ^ ^ + • The scaling 

operators A 0 , a', a" can be used. The subgroup distribution f" 

forms a cluster of many high-order distributions. By random encoun- 
ter, these distributions are homogenized, lose their identity in 
v-dependence and their role of memory-transmission. This random 

behavior leads to a closure of the hierarchy. 

By a macrogroup dynamics, the transport equation for f' is in- 
tegrated to form the "eddy collision" so as to derive the kinetic 
equation for the evolution of the macro-distribution f° in the form 


(3 t ♦ A°L)f° = - L°f ♦ 3*D'-3{f°(t-T)} 
The "collision operator" 


(2.3) 


C'o = 3 • d’*3 (} 


(2.4) 
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depends on the "diffusion operator" d' due to e' -fluctuations. 

By taking the moment, the kinetic equation is transformed into a 
"renormalized Navier-Stokes equation" 

(9 + u«V -vV 2 )u° = E° ♦ J° , (2.5) 

i.e. in a form identical to the original Navier-Stokes equation aug- 
mented by the "eddy damping" 


J° 

u 


dv v c'{f°(t-T)} 


( 2 . 6 ) 


The eddy damping represents the coupling between the macrogroup f° 
and the microgroup fluctuations that organize the collision 
operator. At the same time the equation of spectral evolution is 
obtained as 


h 9 < u 02 > = - T° 

t w. u 

This equation involves the transfer function 


(2.7) 


T° =- <u° • J° > 
u u 


i 2 

K (Vu°) direct cascade 
~x'< u 02 > inverse cascade 


( 2 . 8 ) 


that governs the cascade transfer by direct and inverse cascades 
across the spectrum. By a microgroup dynamics, the eddy viscosity 
k' and the coefficient of eddy amplification are determined from 
the collision operator of the kinetic equation. The change of 
signs indicates the change of eddy damping into eddy amplification 
so that - k'{ J may be interpreted as an operator due to "negative 
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viscosity" . 

t .. t 

The approach of the transport properties K and A to 
equilibrium is regulated by a relaxation process as described by 
the path perturbation £"(t) in a time interval t along the trajec- 
tory that is perturbed by turbulence. The path perturbation enters 
into the orbit function < exp - ik" • i"(T) > and is related to the 
path diffusivity k'^ . 

In order to replace the customary mixing-length hypothesis 

needed in modelings of turbulence (Monin and Obukhov, 1953; Mellor 

and Yamada, 1974, 1982), we formulate a subdynamics of relaxation. 

It consists of a Fokker-Planck equation of path transition and a 

Langevin equation of turbulence. The former determines the orbit 

1 1 

function in terms of the path diffusivity K^,. The latter derives 

li 

the relationship between the path diffusivity K, and the field dif- 

n I 

fusivity D . Since K has been previously derived from the 
D' -dependent collision operator c' in the microgroup dynamics, we 

I! 1 

obtain a system of integral equations for and K . Their solu- 
tions determine the eddy viscosity K* that approaches equilibrium 

ll 

by relaxation from K„. 

The eddy viscosity is derived as 

k' = c K R 0 '* k -2 (2.9) 

With R° = < ( V u° ) 2 >, c., = (2/9)^ P + = 0.471, for P = 0.545 . 

.. — K x z 

The spectral functions of energy, temperature variance, and pres- 
sure variance are found as follows: 

F(k) = Ae 2//3 k -5 / 3 , F (k) = Bc q c- j/3 k" s/3 , F (k) = Ct 2 £ h/i k -7 / 3 (2.10) 

DO P 
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with A = 1.650, B = 0.899 and C = 0.907. The rate of dissipations 
are £ and . 

It is to be noted that a singlet distribution f°, that is 
governed by the kinetic equation of turbulence with the diffusivity- 
dependent collision operator, suffices for the determination of the 
spectrum by our group-kinetic kinetic method, while a system of two 
kinetic equations for the one-point distribution F^t.x.v) and two- 
point distribution F J2 > X 1 » v j ;t 2 » x 2 » v 2 ) is necessary by the 

__ W* 

Bogoliubov method. 

The detailed description of the group-kinetic theory is pre- 
sented in Technical Report #1. 

By the group-kinetic method, the Navier-Stokes equation is trans- 
formed into an equation with "eddy damping" or "eddy am- 
plification". Such an equation is called the renormalized Navier- 
Stokes equation. Since the eddy coefficients are derived 
analytically, the renormalized equation forms a logical basis for 
the modeling of turbulence. Such a modeling is superior to the 
customary modeling by the artificial damping with the differential 
V n of n-th order or by an empirical second-order closure. 

It is to be remarked that the determination of the transfer func- 
tion from the collision operator in the group-kinetic theory in- 
volves the quadruple correlation < E'(t } e'(. t-x) > < f° (t) f° (t-x)> in the 
form of a product of two binary correlations, and not in the form 
of the quadruple correlation < £(t) ^(t-x)u°(t ) u° (t-x) > that can only 
be factorized by the hypothesis of quasi-normality in the perturba- 
tion expansion theory. 
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The group-kinetic decomposition into transport processes yields 
the property of "nearest-neighbor group-interaction". This 
property leads to many simplifications in the development of the 
theory. 


3. Group-kinetic theory of Rossby wave turbulence 

To include the Coriolis force in the Navier-Stokes equation 

A A A 

(2.1), we write E ■ E + E^ as consisting of a pressure field 

E = - — V p and a Coriolis force E c = - X i = - E c + E 8 ’ that 

A A A A A A 

is the sum of E c = u- 0 u x : and Eg= B x 2 u x z . The Coriolis 


parameter oj 0 + 8 x 2 

ferential rotation B . 


has a uniform rotation u. and a dif- 


The curl differentiation ; = V x u of the Navier-Stokes equa- 
tion yields the well-known vorticity equation for a drift wave 


(8^ + u • V - vV 2 K = 


Bu. 


(3.1) 


in two dimensions x = (xj, x 2 ). 


By dimensional arguments based upon the parameter B, the energy 
spectrum 


a 2 l- - 5 


F 2 2 (k J % B" k 


(3.2) 


from u -fluctuations in the drift-range has been proposed by Rhines 
(1975) . 


In the inertial-range that is not controlled by B, Eq. (2.1) is 
reduced to the form 


(a t + u- v- vv 2 )c = o 


(3.3) 
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that identically governs quasi-geostrophic turbulence. Quasi- 
geostrophic turbulence has been treated by Kraichnan (1967), Her- 
ring (1980), and Tchen (1982, 1983), finding the spectrum 

F(k) - CeV* r 1 (3.4) 

It has been tacitly assumed that the uniform rotation u 0 plays 
the role of a circular translation and will not influence the eddy 
transport, so that Rossby wave turbulence without drift becomes 
reduced to the guasi-geostrophic turbulence. However, evidence in- 
dicates that cj 0 controls the eddy transport coefficients. Hence 
the loss of u o -dependence in the vorticity equation (3.1) is a weak- 
ness of the current vortex dynamics because the Navier-Stokes equa- 
tion has been prematurely differentiated. 

It is true that the large-scale structure is drifted by the dif- 
ferential rotation B as indicated by Eg. (3.1), but for the main- 
tenance of the large-scale structure, the small-scale transients 
must organize a ^-dependent eddy amplification in order to perform 
the inverse cascade. Without this transfer, the drift cannot find 
the necessary balance to produce the spectrum (3.2). On the other 
hand, the transients can organize an eddy damping by the 
u) 0 -dependent eddy viscosity to perform the direct cascade in the 
inertial range. Finally the change of the direct cascade into the 
inverse cascade defines a transition range in the spectral 
distribution. The three spectral ranges (drift-range, inertia- 
range and transition-range) cannot be formulated from the 
traditional vorticity equation (3.1). Thus a new vortex dynamics 
is necessary. 
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By the group-kinetic method, we derive the renormalized vor- 
ticity equation 


(3 + u° * V - vV 2 -C ' );° = 

t w. < 


- B U° 4 j° 

2 C 


(3.5a) 


It has all the terms of the traditional vorticity equation (3.1), 
augmented with the eddy damping 


J" = Vx , J° 
£ u ’ u 


dv v C' { f 0 (t -t) } 


(3.5b) 


Upon multiplying Eq. (3.5a) by C and averaging, we find the equa- 
tion of evolution for the enstrophy 


j 3 t <; 02 > = k c - t“ 


(3.6) 


to be governed by the drift function 


K° = 6 2 Kj 

a 


(3.7a) 


and the transfer function 


•p 0 __ 

C ’ 


K'<(V£ 0 ) 2 > direct cascade 


-X '< G ° 2 > 


inverse cascade 


(3.7b) 


The eddy coefficients are found as 


K' 


.5 _ i 

— LO 

2 o 


11 M 

dk F(k ) 


3 

— U) 

4 o 


1 r' , r' = 2^°° dk" k" 2 F(k") 


K°. - 2 u " l 


dk* F,,(k') 


J o 


(3.8) 
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and are proportional to u> 


The eddy viscosity 


k’ 


has been suggested by Blackadar (1962) empirical- 


The energy spectra and the eddy coefficients are found: 


F S ,(VQ = t B* k 


2 v-5 


(3.9a) 


X = oj 0 " 1 8 2 k- 2 


(3.9b) 


(ii) inertial-range 


F(k) = (8/3)'* (u 0 tX k- 3 


K' = (3/2)^ (£_/<-)** k‘ 2 
£> o 


(iii) transition-range 


F(kj = 2(u 0 e ; /8 2 )k- 1 

Here is the rate of enstrophy dissipation. 


(3.11) 


By comparing (3.10a) with (3.4), it is seen that the inertial- 
ranges of Rossby wave turbulence and quasi-geostrophic turbulence 
share the same -3 power law but differ in intensities by the ratio 


K/e ; ) l/6 > 1 

in strong rotation. 


(3.12) 


The detailed theory is presented in Technical Report #2. 
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4. Renormalized vortex dynamics 


By taking the curl differentiation of the Navier-Stokes equation 
we obtain the following equations for the evolution of the vor- 

A 

ticity ; = V xu : 

( i ) vorticity equation for Rossby wave t urbulence — in three 
dimensions 


(3 + u*V - vV 2 ); * (C*V) u - 6 u 2 z +oj 0 V 3 u, V*u =0. (4.1) 

(ii) Vorticity equation for quasi-qeost rophic — turbulence — in_ 
three dimensions 

(3^ + u*V -vV 2 )£ = (C’V)u « ^* u = ® • (4.2) 

These equations have been considered by Lundgren (1982) and Saffman 
(1985) . 

In two dimensions, Eq. (4.1) and (4.2) are reduced to: 

(iii) Rossby wave turbulence 

(3 + C»V - vV 2 ); = - Bu 2 , V-u = 0 (4.3) 

t V W* ■* 

( iv) cruasi-geostrophic turbulence 

- (4.4) 

(3 4 u*V - vV 2 )C = 0 • 

Equations ( 4 . 3 ) and ( 4 . 4 ) are the traditional equations for Rossby 
wave turbulence and quasi-geostrophic turbulence. The losses of 
pressure-gradient I V P and uniform rotation u 0 have been seen as 
expedient simplifications for analysis and numerical computations. 
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The advantage is however illusory. The losses of u 0 and of Vu as 
entailed from Vp have prevented the formulation of the eddy coef- 
ficients that depend on those parameters. 

/s 

With regard to (4.4), the omission of Vp changes the numerical 
coefficient but not the power law k" 3 . But the loss of the eddy 
coefficients that depend on Vu misses the analysis of 
baroclinicity. 

In order to avoid these difficulties, we develop a renormalized 
vortex dynamics and exploit further the group-kinetic method. This 
program will be investigated in future research. 

5. Surface boundary layer 

By the group-kinetic method, we have derived the spectral struc- 
ture in an atmospheric surface layer. In the stable layer, the 
spectrum takes the form k* 5 / 3 , it -1 , gap, k -3 , due to isotropic tur- 
bulence, shear turbulence, stable convection, and Coriolis 
rotation, respectively. In the unstable layer the spectrum takes 
the form k' 5 / 3 , k' 1 and k -3 due to isotropic turbulence, shear tur- 
bulence and Coriolis rotation. The gap disappears in unstable con- 
vection. These results are verified by experiments. The 
detailed theory is presented in Paper #3. 


12 


6. Two-phase turbulence 


Consider a particulate phase of number-density N velocity u 2 , 
pressure p 2 suspended in a fluid phase of velocity Uj pressure Pj 
and constant density p . The two phases are coupled by the constant 
friction y. The governing equations are: 

N _1 p£(3 t ♦ Uj ’VJUj + ^ VpjD = y(u 2 - Uj) 

V • u = 0 

•v ** * 

m s C3 t 4 u 2 *V)u 2 + Vp 2 =-y(u 2 - Uj) (6.3) 

^ M V A ^ 

3N*V*(Nu)«0, ( 6 • ^ ) 

^ w v. 
r\ 

For p 2 a constitutive equation is added. 

The two species are denoted by ( ) with a = i, 2 for fluid phase 
and particulate phase, respectively. An instantaneously fluc- 
tuating function is denoted by ()=()+() and can be decomposed 
into an ensemble average (")=<> and a fluctuation (^). 

The variable number density N does not enter in the momentum 
equation for the particulate phase, but does enter in the momentum 
equation (6.1) for the fluid phase. This equation resembles the 
one that governs convective turbulence, where the Boussinesq approx- 
imation is acceptable if the phase-coupling can be compared with 
the convection. With this approximation the number density in the 
left side can be assumed to be a constant N , simplifying the 
system of equations (6.1) - (6.4) into 


( 6 . 1 ) 

( 6 . 2 ) 
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(6.5) 


K 0 ' 1 EPO t - Uj-VjG, + V Pi 3 = Y (u 2 - u t ) 

V*Uj=0 (5.6) 


m c (9 t + V)u + V P, s - Y(u - u ) 
S I 2 2 ^ * l 

V • u 2 = 0 


(6.7) 

( 6 . 8 ) 


If VP << v(u 2 - Uj), the system of equations (6.1) - (6.4) is 

reduced to 


N- , [P(3 t + u j *V)Uj ♦ pj] = Y(u.- Uj) (6.9) 

V • Gj = 0 (6.10) 

m s (9 t + u-V)u 2 = - y(u 2 - u,) (6.11) 

3 t N «■ V.(N G 2 ) = 0 (6.12) 


for dilute suspension, called "dusty gas" by Saffman (1962). 

The treatment of two-phase turbulence from the system of equa- 
tions (6.5) - (6.8) has been our first attempt. The spectral re- 
sults are: 
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a. In the inertia-range 

F (k) = A c 2 / 3 k" s/3 (6.13) 

a a a x ' 

with c a = v a <(Vu a ) 2 > + a a <\ 2 >* = n 0 Y/p ' q 2 * Y/m s » vhere ( ) b i s 

the other species. 


b. In the drag-range 


F ‘ (k) = B i U I X k “ > » F (k) - B u 2 k- 1 

^ 2 *2 


(6.14) 


with u = (e /a )' i . 

*a v a a' 


The numerical coefficients a , B are 

a a 


evaluated. The detailed description is presented in Tech. Rep. #7, 


The treatment of the dusty gas model (6.9) - (6.12), as well as 
the compressible model (6.1) - (6.4) will be done in future re- 
search. 


7. Soliton turbulence 


The soliton turbulence is governed by the Schrodinger equation 

(i 8 ♦ v V 2 - x u A E 2 )E = X . 
t n 2 n * 1 _ 


Here S' is the envelope of field fluctuations, such that the density 

0/ 'V, 'V 'V 

is N % A | E | , X is the driving force, and are constants, and 
A = 1 - A is a fluctuating operator and gives the deviation from the 
average A= < >. This equation applies to nonlinear problems in hy- 
drodynamics (gravity waves and surface waves, shallow water waves) 

optics (nonlinear propagation, self-focusing, saturation of scintil- 

. . 'v 

lation) and biological systems. Thus E depends on the nonlinear 
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problem considered. The transformation of the hydrodynamic equa- 
tions by taking the envelope of the wave packets is made (Tchen, 
1986a). We derive the spectral structure F^ ^ k -1 and F^ k 5 f° r 
density fluctuations, and 

F * k" 1 , F c ^ k -3 

for field fluctuations. The results agree with plasma experiments 
(True, 1984). The detailed description is presented in 
papers #8, #9. 

8. Compressible turbulence 


The following models may be used for the description of the 
microdynamic state of turbulence. 

(i) The Kavier-Stokes equations are 


3 P u. + V.pu.u. « pE . , with pE = - Vp 
t l ) j l pi P 


/V ''A 

3^p + V*pu = 0 


The equation of momentum may be written as 


(8.1a) 

(8.1b) 


(3 + u • V)u = E 

t _ - - -P 

A A A 

A constitutive relation p = pRT for ideal gas or 

A A A 

adiabatic gas can be used. Here p, u , E, R , y are 
sity, velocity, pressure gradient, gas constant 
specific heats, respectively. 


( 8 . 2 ) 

(p/P 6 ) = Cp/p 0 ) Y for 

pressure, den- 
and ratio of 


(ii) The equations of sound propagation in turbulence are 
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(3 + u • V)u = E 

(3 2_ V-c 2 V) p = r 

L v 


(8.3a) 

(8.3b) 


with a source r = V V: puu . The equation of propagation (8.3b) 

As ^ ^ y 

with variable speed c = (dp/dp )” 2 is obtained by a cross differen- 

tiation of (8.1a) and (8.1b) with respect to and t. 

(iii) For an adiabatic gas with 



and by introducing the speed of sound 
c = (Y p /PK 

and an auxiliary function 



We transform Eqs. (8.1a) and (8.1b) into the following system of 
symmetric equations 


(3 t + u • V)u 


111 

2 


K • V K 


(8.4a) 


~ ~ v - 1 ~ 

(3^ + u • V)K = K V *u 


(8.4b) 


This system resembles the Riemann equations for rarefaction 
waves. Note that 




.2 


Y-l 


gives the potential of the driving force 


(8.5) 
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( 8 . 6 ) 


E = - V4> 
_P - P 


(iv) For interstellar gas, the hydrodynamic equations are 


^ ^ AA A A 

9..PU. + V.pu.u. = E. 
* i 111 i 


A A A 

3^p + V • pu * 0 


The driving force 


E = 


E 

-P 


+ E 
~g 


(8.7a) 

(8.7b) 


( 8 . 8 ) 


consists of a pressure gradient 



and a gravitational force 

A 

E - - 
.g ~ g 

A 

where the potential <t satisfies the Poisson equation 


(8.9a) 


(8.9b) 


V 2 4> g = 4 ttGp (8.9c) 

For the problems of compressible turbulence above, the group- 
kinetic method is best suitable. By raising to higher dimen- 
sionality by 

A ~ 

f(t,x,v) = p 6 [ v -u(t ,x) 2 , (8.10) 


Eqs. (8.1a) and (8.1b) is transformed into the master equation 


(3 t + v *V + E • c) f(t ,x,v) = 0 


(8.11) 
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It is not difficult to show by moments J dv f = p , J ^ * = p Jf. 

that the master equation can reproduce the Navier-Stokes 

equations. An instantaneously fluctuating function as denoted by 

~ % _ 

(J = t) + (.) can be decomposed into an ensemble average ( ) = < > 

and a fluctuation ( ). It is interesting to see that turbulence, 

whether incompressible or compressible, is represented by the same 

master equation for the description of the microdynamic state of 

turbulence. 

We have investigated the Navier-Stokes equations (8.1a) and 
(8.1b) for compressible turbulence by the group-kinetic method. 
The following spectral intensities are derived: 

< u' 2 > % k' 2 , < e' 2 > ^ k" 2 , < p' 2 > % k' 1 * . (8.12) 

The details of derivation are given in Technical Report #13. 

9. Conclusions 

Past analytic efforts in the literature have not succeeded in 
deriving explicit expressions of cascades (direct and inverse 
cascades) and of eddy coefficients (eddy viscosity, eddy coef- 
ficient of amplification, and eddy diffusivity of drift), that con- 
trol the spectral evolution of turbulence. Since the spectral 
evolution, the transport coefficients and the relaxation for the 
approach of the eddy coefficients to equilibrium are the three tran- 
sport processes of turbulence, we represent them by a macro-kinetic 
group, a micro-kinetic group and a subgroup, respectively. 

By group-kinetic considerations, we formulate a kinetic equation 
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of turbulence, and derive the eddy coefficients analytically. The 
approach to equilibrium is analysed by a subdynamics that consists 
of a Fokker Plank equation of path perturbations and a system of 
Langevin equations of turbulence. We derive the "renormalized 
Navier-Stokes equations" and the "renormalized vorticity 
equations". They contain explicit cascades through the 
analytically derived eddy coefficients. 

By the group-kinetic method we have investigated some aspects of 
Navier-Stokes turbulence with the Kolmogoroff spectrum, Rossby wave 
turbulence, boundary layer turbulence, two-phase turbulence, 
soliton turbulence, and compressible turbulence. Other aspects 
will be investigated in future research. 
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GROUP-KINETIC THEORY OF ROSSBY WAVE TURBULENCE 


By C. M. TCHEN 

Department of Mechanical Engineering, City College of New York 

New York, NY 10031. 


The Rossby wave turbulence is formed by large-scale two-dimensional 
vortices in the atmosphere or ocean on a rotating planet. It has 
long been thought that the standard vortex model 
Q +u*V-vV 2 )c = -6u is the basic governing equation. Here u is 

t - — 2 

fluid velocity, z, = Vxu is vorticity, v is molecular viscosity, and 
£ is dif ferential, rotation. The present wor)t demon&strates that 
this model would fail in analysing the spectral composition of the 
Rossby wave turbulence, because it had canceled the uniform rota- 
tion u 0 upon which the eddy transport coefficients must depend. A 
renormalized vortex dynamics is developed by the use of the group- 
kinetic method. It describes the evolution of the large-scale vor- 
tices under drift and, in addition, their w o -dependent interactions 
with small-scale transients. The velocity distribution is decom- 
posed into three groups representing three transport processes 
(evolution, eddy transport coefficients and relaxation) . The 
kinetic equation of turbulence is derived. The eddy coefficients 
are calculated from the "collision operator". A subdynamics deter- 
mines the relaxation to equilibrium, and consists of a Fokker- 
Planck equation of transition and a system of two Langevin 
equations. The closure is found by memory-loss and not by the 
hypothesis of quasi-normality. 
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The spectral structure is divided into an inertial-range by 
direct cascade, a transition-range from direct cascade to inverse 
cascade, and a drift-range by inverse cascade, in the order of 
decreasing wavenumbers. The corresponding spectral laws are found 
to be: F(k) = Aj (w # e^k' 3 , F(k) = e /B 2 )k‘ l , F 22 (k) = A^k” 5 . 

L 

The numerical coefficients are evaluated: Aj = C S/3) , A 2 * 2 and 

A 3 = 8/3. 


1. Introduction 

The Rossby wave turbulence is described by the Navier-Stokes 
equation of motion with a Coriolis force that has a uniform rota- 
tion io and a drift from the differential rotation 6. The works on 
0 

large-scale turbulence are well documented (Charney, 1971; Baer 
1972; Salmon, Holloway and Hendershott, 1976; Salmon, 1978; Lam- 
bert, 1981) and have shown a spectrum k“* in the inertial-range. The 
change to the spectrum k -5 in the drift range was predicted by 
Rhines (1975). Observations have shown a spectrum k 1 in the tran- 
sition range. It will be desirable to distinguish between the 
Rossby wave turbulence and the two-dimensional geostrophic tur- 
bulence in their inertial-ranges where the same power law -3 is 
shared. 

For the description of the Rossby waves that are large-scale 
planetary vortices in the rotating atmosphere and ocean, a vortex 
dynamic is required. The curl differentiation of the Navier-Stokes 
equation 
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( 1 . 1 ) 


( 3 ♦ u • V - vV 2 ) u = E , 


with 


^ A 


E 4 E , E 
-P -c -p 




A ^ ^ 


U 0 UXZ, Eg e Bx 2 UX Z, 


“ c e “V 6x 2 * 


** ( 0 , 0 . 1 ), 


yields the vortex model 


(1.2a) 

(1.2b) 


O t ♦ u«V-vV 2 ); * -$ u 2 (1-3) 

for the evolution of the vorticity £ » Vxu in two dimensions with 
a drift by B. The fluid velocity satisfies the condition of incom- 

A 

pressibility V*u = 0. Here v is molecular viscosity, p is pressure, 
p is density. The^ Coriolis parameter w has a uniform rotation u> 0 
and a differential rotation B. A fluctuating function, as denoted 
by () c ()+() is the superposition of the ensemble average 
() = < > and the fluctuation C) • 


In the inertial-range of the Rossby wave turbulence, the drift 
is negligible, so that Eq. (1.3) is reduced to the homogeneous form 

(9 t ♦ u • V -v V 2 ); * 0 (1.4) 

that is identical to the equation for two-dimensional quasi- 
geostrophic turbulence. This simplification has appealed to theore- 
tical studies and observational interpretations, but is deceptive, 
especially since u 0 should control the eddy coefficients and 
governs the spectral structure of the Rossby wave turbulence. 

We believe that a proper vortex dynamics should describe the 
evolution of the large-scale structure by the curl differentiation 
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of the renormalized Navier-Stokes equation that includes the in- 
teraction with the small-scale transients for organizing the 
io 0 -dependent eddy damping. In this respect, recall that a 
u> 0 -dependent eddy viscosity has been suggested by Blackadar (1962) 
empirically. 


A particularly appealing method for treating turbulence has been 
the "renormalization perturbation expansion". It has enabled 
Kraichnan (1959, 1977) to develop his DIA method and 
generalizations. A systematic method of developing the expansion 
was given in a general form by Martin, Siggia and Rose (1972) and 
by Dubois and Espedal (1978). Without a convergent summation, the 
method of perturbation expansion cannot easily derive the eddy coef- 
ficients which may be of kinetic origin. The couplings among modes 
in expansion are too complicate and cannot clearly describe the 
physical processes they represent. For this reason, we develop a 
group-kinetic method for Rossby wave turbulence, not from the 
vortex model (1.3) but from a "master equation". The master equa- 
tion is based upon the Navier-Stokes equation and describes the 
microdynamic state of turbulence. It is only after having deter- 
mined the u 0 -dependent eddy transport coefficients and derived the 
"renormalized Navier-Stokes equation", i . e . a Navier-Stokes equa- 
tion with "augmented damping", that we pass to the "renormalized 
vorticity equation". Our vortex dynamics will have the same form 
as Eq. (1.3) but is augmented with an "eddy damping" that depends 
on u° and is analysed from the "eddy collision" of the kinetic equa- 
tion of turbulence. 


By group-kinetic scaling, the three transport processes of spec- 
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tral evolution, eddy viscosity and relaxation are represented by 
three groups of distributions (Section 2). In the macrogroup 
dynamics, we transform the Navier-Stokes equation into a master 
equation and derive the kinetic equation of turbulence (Section 
3). The "diffusion operator" with orbit functions enters into the 
"collision operator" (Section 4). The kinetic equation is reverted 
to the renormalized Navier-Stokes equation by moment, and subse- 
quently to the renormalized vorticity equation by curl differen- 
tiation. The evolution of enstrophy and the governing transport 
functions are found (Section 5) . 

In particular we are interested in a microgroup dynamics ' .at 
derives the eddy transport coefficients from the collision oper. ror 
(Section 6) . 

A subdynamics of relaxation is developed in Section 6 for ~he 
approach of the eddy viscosity to equilibrium, and consists o: a 
Fokker-Planck equation of path transition and a system of rwo 
Langevin equations of turbulence. The subdynamics determines the 
orbit functions and the path diffusivity. Finally, the transport 
coefficients (eddy viscosity, coefficient of amplification, and 
drift diffusivity) for the direct cascade, the inverse cascade, and 
the drift are investigated in Section 7. The spectral structure is 
found in Section 8. A summary with discussions is presented in Sec- 
tion 9. 
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ORIGINAL PACE IS 
Of POOR QUALITY 



2. Group-kinetic method 


Being a nonhomogeneous partial differential equation, the Navier- 
Stokes equation is not convenient for the formulation of a theory 
of transport in the gradient form or in the form of an eddy damping 
(or amplification) . We raise its dimensionality and transform it 
into a homogeneous equation 


O t ♦ L)f = 0, 

called the "master equation" 


( 2 . 1 ) 


by writing f(t,x,v)= 6Cv -u(t,x)D . Here L=vV-vV z «-E*3, 3 t =3/3t, 3=3/9v. 
Since v is an independent variable the master equation has lesser 
nonlinearity. It is # not difficult to verify by taking the moments 
that the master equation reproduces the Navier-Stokes equation and 
the equation of continuity. 


From the statistical viewpoint, the master equation describes 
the microdynamical state of turbulence where the elementary interac- 
tion among fluid elements is represented by £ -fluctuations. A 
kinetic theory can be developed, either by the reduction of the 
N-point distributions into F Jf F 2 , F s according to the method of 
Bogoliubov (1962), or by the decompositions 
%f°+f ' , f ' f"=f + - • • from the group-kinetic method (Tchen, 
1978, 1984, 1986). The three groups f 0 , f * » f" called macrogroup, 
microgroup and subgroup represent the transport processes of spec- 
tral evolution, eddy viscosity, and relaxation, respectively, with 
durations of correlation t° > t' > t" in the order of increasing ran- 
domness. The scaling operators A 0 , A*, A" can be used. 
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We choose the group-kinetic method in view of its benefits over 
Bogoliubov's method and the method of perturbation expansion. For 
the derivation of spectrum, the singlet distribution f° suffices in 
our group method, while two distributions F , F 2 are needed in 
the Bogoliubov method. The eddy viscosity being of kinetic origin 
cannot be easily derived by a continuum theory. Although the 
Fourier modes in groups interact indefinitely and may overlap in 
wavenumbers, their statistical properties are well separated by the 
transport processes they represent. Consequently, the groups in- 
teract statistically as nearest-neighbours . Hence a quadruple cor- 
relation takes the form of a product of two binary correlations as 

< f°(t)f°(t-T) > < E'(t) E'(t-T) > (2.2) 

V 

and not in the form < u(t)u(t-T)l’(t)^(t-T) > that requires the 
hypothesis of normality for factorization, as was in perturbation 
expansion theory. The subgroup f" = . . forms a cluster of 

many high-order distributions. By random encounters, they are 
"homogenenized" and lose their identity in v-dependence. As a re- 
sult, they can be simulated by an effective medium of friction coef- 
ficient C"(k) . Since it is by the v-dependence that the memory is 
transmitted, the loss of v-dependence causes the loss of memory, 
leading to irreversibility and closure, and providing a mechanism 
of relaxation for the approach of the eddy viscosity to 
equilibrium. 

The correlation functions of macrogroup, microgroup and subgroup 
have the limits (0,k), (k,°°) and (k" ,°°) in integrations 

dk’"... • These limits may be regarded as the 
statistical demarcations of the groups. 
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3. Kinetic equations of turbulence 


By the use of A® and A 1 we transform the master Eq. (2.1) into 
the following transport equations 

O t *A # L)f # «- L°f ♦ C 1 { f 0 ) » with C l (f°) = - A® L’f 1 
(9 t *A , L- C'^f 1 «= - l^f® . 

The auxiliary operator A 1 gives A J f = f 1 . Note that f 1 
a medium that offers a friction C’^ 1 . The integration of Eq. 
(3.2) gives f 1 and a multiplication by L 1 yields the eddy collision 
in the form 


(3.1) 

(3.2) 

evolves in 


CMf 0 ) = - A 0 ! 1 f 1 = 3* D 1 *3 { f°(t-T) } (3.3) 

or 

= - A® l' f ' = 3 • D 1 • 3 { f°(t-T) ). (3.4) 

The approximation C 1 =C is a consequence of the nearest-neighbour 
group-interaction. The diffusion operator is 


ft“® 00 

D'{ } = I dT<E'(t) E ' (t -T) >{ ) ^ 

* Jo - (3.5) 

and the collision operator is C ’{ )= 3 • d' • 3 {) . Hence we derive 
the kinetic equation of turbulence 


(3 t +A°L)f° = - L°f + c'{f Q ) . 


C3.6) 
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4. Orbit functions in strong rotation 

The diffusion operator (3.5) is obtained by the time integration 
of the Lagrangian correlation of field-fluctuations. In view of the 
separation of transport processes, the diffusion operator is 
adiabatic, i.e. t-+«°. In strong rotation, i.e. E = E £ , the Lagran- 
gian correlation in the operator description 

< E^(t)E^(t-T) > - < E^(t)A’U(t,t-T)E^(t-t) > , (4.1a) 

or, equivalently, in the orbital description 

< > * < E^(t)E;tt -T, x ( t- T )3 > f (4.1b) 

can be calculated from the evolution operator U(t,t-t) that is 
governed by the equation 

O t ^A , L l -C ,, )U(t,t’) « 0, with U(t,t) • 1. t’<t, (4.2) 

or from the perturbed trajectory 

X(t-T) « X - dT V (T ) - f.(T) 

— — Jo — c — 

* x - mv* z - £(t) (4.3) 

with 


v c (t) k h ti ) 0 T V X z , m=!a w 0 T . 

We choose the orbital description. 

By Fourier transformation, the Lagrangian correlations are 


(4.4) 
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< u°(t)f°(t-T) > * fdk’ X< u e (k') f° V)> h (T,k\v)h £ (T,k') 

(4.5a) 

< E;(t)E^(t-T) > • Jdk" X< E^(k")E^(-V") > h c ('t.k , .v)l> 1 (t.k") 

(4.5b) 

< E”(t)E"(t-T) > - U’"J<< E^k" V>E'") ’) . 

(4.5c) 

the diffusivities are 


D’(t,x,v) - j°°dT jdk >!<E c (k ) E c ( -k ) >h c (l,k ,v)h^(l,k ) 

(4.6a) 

D "(t.x) = j'dT fdk" tr^<E;(k , ")H , ( ;(-k" , )> h^T.k"’). 

(4.6b) 

and the orbit functions are 


h (T,k",v) * expC-imk"* (vxi) J 
C - ^ *•> ^ - 

(4.7a) 

h t (T,k") = Cexp c’ (k ")t3 < exp-ik"*£"(T)>, h'^(T ,k")=<exp-ik* (t)> (4.7b) 

for streaming and relaxation. The path fluctuation £ 
path diffusivity 

organizes a 


* ’s lim^- tr <£"(T)i"(T) > 

= f d-r’ tr < u'(0)u’(t)> . (4.8) 

J o *" 

The meandering caused by £° and S. 1 is the lowest-order perturbation 

«s 

expansion in quasilinear theory and is negliglible in strong tur- 
bulence. 
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The Fourier decomposition is truncated in a region within which 
the function is quasi-homogeneous, and the truncation factor is ^ . 
The limits of the volume integrations in (4.5a) - (4.5b) and m the 
following are understood to extend from - » to + 

5. Renormalized equations, transport functions and eddy 
coefficients 

5.1. Renormalized equations 

By taking the moment of the kinetic equation (3.6) we get the renor- 
malized Navier-Stokes equation 

(9 u° • V - v V 2 ) u® = E° ♦ J* (5.1) 

with the "augmented damping" or "eddy dampings" 

J® = J dv vc'{f°(t-T)} . (5.2) 

By a curl differentiation £® = VXu® we transform Eq. (5.1) into the 
renormalized vorticity ecjuation in two dimensions 

(9 + u° • V - v V 2 ) C° = - & u° * J® . with J® = V X J® . (5.3) 

The eddy dampings represent the coupling between the macro- 
distribution f® and the collision operator C that depends on VP 
and iDju'xz. Thus the pressure gradient and the uniform rotation 
which were missing in the traditional vorticity model (1.3) are 
recovered in the renormalized vorticity equation (5.3) . 

Upon multiplying Eg. (5.3) by &® and averaging, we find the equa- 
tion of evolution 
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(5.4) 


h 3 < ;® 2 > = - B < u® ;•> - v < (V ;• ) 2 > -t® 

t 2 — S 

for the enstrophy *:<£ 02 >. The transport functions are: the drift 

function 

W° = - 6 < uj ;° > (5.5) 

the transfer function 

T® = - < ;° J®>= - < ;°v *J® > (5.6) 

and the dissipation function v<(Vc®) 2 > . 

5.2. Drift function 

The flux -<u°£°>i 5 calculated . by integrating Eq. (5.3) for the 
evolution of 5° rewritten as 

O t 4u° • V-vV 2 - c ’)c°« -BuJ (5.7) 

with a friction constant c'(t,x) . An integration with respect to 
time along the perturbed trajectory, a multiplication by u° and an 
ensemble average yield the flux 

-<u°;°> = 6 Kj , (5.8) 

transforming (5.5) into the drift function 

W° * B 2 K® . (5.9) 

a 

The drift diffusivity 
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K d “ dT <uj (t) U 2 ° (t-T) > 

J Q 

• j dT | d i' X <u°(k')u°(-k’) > h £ (T,k') . (5.10) 

is controlled by the orbit function 

hj^T.k') = C exp c'(k') 3 < exp -ik’ (t) > (5.11) 

for relaxation. 


5.3. Transfer function 

The transfer function (5.6) can be written in the form 
Tj ■ - V *jdv V (3-„D’*3){ <;°(t)f°(t-T) >} 


jdv v 


dk" dk , X<E’(k")E;(-k") > h^T.k^h^fT.k’) 


dx M°(T,k”,k',v){ X <; 0 (t,k')f 0 (t,-k\v) >) 


(5.12) 


by (5.2), (3.4), (4.6a) and by Fourier transformation. The memory 

operator is 

M°{) = 3 - h (t ,k", v) 3 .h (T ,k', v) {} . (5.13) 

the differentiation 3^ 3^ there are terms which cancel by in- 
tegrations with respect to dvdk u dk' . The remaining terms give the 
following memory operators: 

(i) for direct cascade 


M O = h c (r,k , v)3 2 h c (T,k , v) {} , with 3 2 h (t , k ' ,v) = -m 2 k ' 2 h (T.k'v) (5.14a) 

^ “ -v ■<— C » — « 
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lira M°{) = - tt ? k ,2 > in = Jjco 0 t 2 , 
y=0 


(5.14b) 


by (4.7a); 

(ii) for inverse cascade 

M° { } = C 9 2 h (T,k",v)3h (T.k'.v) (5.15a) 

lim M°{) = m 2 k" 2 . (5.15b) 

v=0 

With the loss of memory by omitting the v-dependence, we tran- 
form the memory operator into (5.14b) and (5.15b), and reduce the 
transfer function (5.12) into 


y 0 _ 


J. dI J 


x m 


dk" dk' trX<E^(k") (-k") > h £ (t ,k") 


r_v’21 


^“ tt2 jX < C (k ) C°(-k’) > 


dk' D’{ - 2 


' k 2 |X<C°(k , )C e (-k') >) 
k " 2 


ID 


D' {m 2 } R® , with R° = < (V ;° ) 2 > 


D’{m 2 k" 2 } <; 02 > , with <; 02 > 


k'r® , direct cascade 


-X'<!; C2 > , inverse cascade . 


jdk’ X<;° (k*) ;°(-k') > 


(5.16a) 


(5.16b) 


Here 


D' . |Jdt fdl" ,rX<E' (k ", E'(-k")> h,( T ,k") 


(5.17a) 
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(5.17b) 


K' = D’{m 2 } = fdk" tr* <E'(k") E 1 (-k") > fdT m 2 (i)h f (t ,k M ) 

J ^ ^ J o 

x’ = K’{k" 2 } * |dk" k" 2 tr 3( <E^(k")E^(-k")> j^dT m 2 (T)h £ (T ,k") 


(5.17c) 


are called path diffusivity, eddy viscosity, and coefficient of am- 
plification, respectively. 


The results (3.4), (5.2) and (5.16a) indicate an equivalence 
C'(t,x,v) = c'(t=»,x) = KV (5.18a) 

between the integral operator C* (t.x.v) { } and the differential 

operator C (t= °°,x) , provided the eddy viscosity is computed from 

the eddy collision of the kinetic equation. In an analogous way, we 
have 


C"(x)=K"V 2 and C"(k)=-k 2 K" . (5.18b 
6. Subdynamics of relaxation 

The principal aim of the subdynamics is to determine the path 
diffusivity k'^ and the orbit function h^(T,k") to find a relaxa- 
tion for the eddy coefficients to approach equilibrium. To this 
end, we consider the Fokker-Planck equation of transition and a 
system of two Langevin equations of turbulence. 
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6.1. Fokker-Planck equation of path transition 


The probability for the path to be in the interval 1 and £*d£ in a 
time interval t is p(t,x;t-t, x-£)d£, or p(T,£)d£ by abbreviation. 
The probability density p(t,£) 
malization and is governed 
1944) . 

3 , ii 3 ^ 

jr P(T.£) -= K £ — - pU,£) 

“ 3£ 2 ~ 

The path diffusivity k" as defined by (4.8) is adiabatic, i.e. in- 
dependent of t and £. The Fourier form is 

§7P( T . k ") * - PlT.k") (6.2) 

The solution 

P(T, k ") •= exp(-k" 2 K"T) (6.3) 

~ (2ti) 

in d ■= 2 dimensions determines the orbit function 

h"(T < k 1 ') = < exp-ik"*£(i) > = [d£ e' 1 * *«p(t.£) *= (2Ti) d p(T ,k") 

e (2u) d p(T,k") = exp(-k ,,2 K^T) (6.4) 

6.2. Langevin equations of turbulence 

The orbit function h £ (T,k") depends on the path diffusivity k'^ and 
governs the relaxation for the approach of the u 0 -dependent eddy 
viscosity to equilibrium. Note that B governs the drift in the 

macrogroup transport (5.9) and does not enter into the subdynamics. 
We consider a system of two Langevin equations 




satisfies the condition of nor- 
by the Fokker-Planck equation (Tchen, 


( 6 . 1 ) 
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du G 


U- * 

g" 

dt 

G 

_G 


du" 

dt 

y G 

«« 

u * 

So- 


The first Langevin equation has 
correlation 


(6.5a) 

(6.5b) 

a random noise g" with an auto- 


tr < G (0 )G"(t) > = 


( 6 . 6 ) 


that depends on the field diffusivity D • It determines the fric- 
tion constant Y q for the second Langevin equation to analyse k£ 
and calculate h^. In the second Langevin equation, the field is 



u„ 


n A 

U G X 2 


(6.7) 


by (1.2b). The system is immersed in a quasi-geostrophic tur- 
bulence, labeled by the suffix ( as the background medium with a 
spectrum and an intensity of the form 


F(k) 


C G e 


and < u” 2 > * C G e^ 2 /*k 


( 6 . 8 ) 


Here * v<(7i;) 3 > is the rate of dissipation and is a numerical 
constant. 
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6.3. Determination of the friction constant Y r from the first Lanqevin 
equation b 

Upon multiplying the first Langevin equation (6.5a) throughout by 

II 

the path length X G and averaging, we have 

d£" 


* (if) * ^ < ys > * 0 


or 


- <u g 2> + 2 ^g k L'° * 


and obtain the relationship between the path diffusivity K and 

xG 


the friction constant y g in the form 


II 

'LG 


* b Y / 1 < u ' 2 > 


' b C G c ; 2/3 V lk ' 2 


(6.9) 


by (6.8). We have written k",. * tr<u r. *r. > 

.2 , .*'2 . 2 


The term 


XG _ G .6 

h d < £ g ‘ > /dt z is negligible in statistically steady tur- 
bulence. 


On the other hand, we calculate 


*= j^dt Jdk" tr X < UgCOugC-k') > h^CT.k”) 

■ Jdk tr X < u G (k")Ug(-k") > J°°dT exp(-k' 2 K^-t) 

= f dk" tr X < UG (k ,, )u G (-k") > (k" 2 ^)- 1 , (6.10a) 

by (4.8) and (4.7b). The integral equation gives the solution 
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II 

K £C 


C 2 



dk 



F(k") 



* Cb C c t c 2 / *) h k - 2 


(6.10b) 


by (6.8). 

By a comparison between (6.9) and (6.10b) the friction constant 
is determined as 

y G * (! * C G c G 2/,)li * (6 ' 11) 


6.4. Determination of the relationship between k" and D" from the 
second Lanqevin equation *■ 

From the second Langevin equation (6.5b), we derive the relation- 
ship 

K £ c V ?D " (6.12) 

between the path diffusivity and the field diffusivity 
D" = | dk" tr X < E’^k") E^-k") > £ dT h £ (T,k") 


J dk" F c (k") 


• ( k " 2 K ")- 1 


by (6.4). The detail of calculations has been omitted, 
trum F c (k)| G has been introduced such that 

[ dk" tr X < E^OE^C-k") > = f dk" F (k")| . 

J - - - J k ' 


(6.13) 
The spec- 


The system of equations (6.12) and (6.13) are solved, determining 


D * [2 y g z | dk k”' 2 F c (k") 
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(6.14a) 



(6.14b) 


K’’ = Y G - l C2 J" dk" k"- 2 F c (k”) 3 


By writing F c (k)| G * u o F ( k ) 
(6.14) and find 


D £ y g 2 u> 0 k- 2 


from (6.7) and (6.8), we calculate 


(6.15a) 


*i = “o k 


-2 


(6.15b) 


It follows 


\ * exp(-w 0 T) 


di h^(T) = w 0 "* f 


o di T* h^(T) = 24 u) o - 5 , 


(6.16a) 

(6.16b) 


by (6.4) 


7. Eddy coefficients 

The eddy viscosity K* and the eddy coefficient X* from (5.17b) 
and (5.17c), determine the transfer function in direct cascade 
(5.16a) and inverse cascade (5.16b), while the eddy diffusivity K® 
in (5.10) determines the drift function (5.9). These eddy coef- 
ficients approach their equilibrium by relaxation through the orbit 
functions h^(T,k") and h^(i,k") which differ by c"(k") = -k” 2 K", 

from (5.11) and (5.18b). However, an inspection of (5.17b) and 
(5.14b) reveals that K" is inversely proportional to while is 

proportional to w 0 by (6.15b), so that the contribution of c" 
becomes negligible, approximating 

hi*"! f 7 - 1 ' 
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7.1. Field diffusivity 


The field diffusivity D (t,x) is defined by (4.6). By writing the 
field intensity trX < E^(k') E^-kV = u 2 tr <u' (k ')u' (-k')> and by the use 
of (6.16b) we get 

d' (t,x) {} = u o jdk" IT x< u , (k , ’)u , (-k")>{) 

= «. dk" F(k") {} . (7.2) 

7.2. Eddy viscosity K* 

The eddy viscosity as defined by (5.17b) is K 1 = D* (t ,x) {m 2 ), 
with m 2 = (’4 U D t 2 X 2 , by (5.14b), or 

k' = (h - fl ) 2 d' {t 1 * ) 

c ( ! < w e ) 2 j^ dk" F(k ,, )j o dT T-h^a.k*') 

e J w o' 1 £ dk " F(k '^ » (7.3) 

by (7.2) and (6.16b). 


7.3. Coefficient of amplification X' 

The eddy coefficient X 1 as defined by (5.17c) is calculated from 
the integral operator K * { ). Thus we find 


X' = j* dk" k " 2 F(k") = |u ( .* 1 R'(k), 


(7.4) 
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.iiJiiii 


from (7.3), with the vorticity function R* = 2 J d k k 2 F(k' ) 


7.4. Drift diffusivity k p 


The drift diffusivity K° is defined by (5.10) and can be cal- 
culated in the sane manner as for (7.3). We find 


K d * 2 


f* 

-o ' 1 f 

J 0 


dk Fj 2 (k ) , 


(7.5) 


f k 

where F 22 (k) is the spectrum such that < u°u° > = 2 dk’ F 22 (k') * 

J n 


The expressions (7.3) - (7.5) show that the eddy coefficients k', 

l’ and K° are proportional to u.' o " 1 . 


i j 


8. Spectral structure of Rossby wave turbulence 

We consider three ranges of Rossby wave turbulence: (i) inertial- 
range by direct cascade, (ii) drift-range by inverse cascade, and 
(iii) transition-range from direct cascade to inverse cascade in 
the direction of decreasing k. 


8.1. Inertial-range by direct cascade 


In the inertial-range, the transfer occurs at the constant rate of 
enstrophy dissipation and is governed by 


K 


by (5.16a), where k’ is given by (7.3) and r® 
In terms of F(.k),Eq. (8.1) is 


( 8 . 1 ) 

k 

dk ' k ' ** F(k’). 
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( 8 . 2 ) 


3 

— Cl) 

2 o 


■> f dk" F(k") [V k’ 1 * 

Jk JO 


F(k') 


Upon dividing both members by any one of the integrals and differen- 
tiating, ve solve the integral equation and find the spectrum 

F(k) - (8/3) ** (oj 0 k-’ (8.3) 

and the intensity 

< u' 2 > = (8/3) 1 * (u 0 k" 2 . (8.4) 

The eddy viscosity is 


K* »= (3/2) ls Cc c /uj o ) 1s k‘ 2 , (8.5) 

by (7.3). The ’’spectrum (8.3) should not be confounded with the 
spectrum of quasi-geostrophic turbulence (6.8) although the two 
spectra share the same power law. 


8.2. Drift-range 

The large-scale structure formulates a drift K° , while the small- 
scale transients perform an inverse cascade to amplify the en- 
strophy at the rate X . The maintenance of the large-scale struc- 
ture is described by the balance between the drift and the cascade, 
a phenomenon called "blocking". The balance in the differential 
form is 

o o 

dK dT 

2- «= 0 . 

* d (" k ) ( 8 . 6 ) 

Here d/dk and d/d(-k) are differentiations in the direction of in- 
creasing and decreasing k, respectively. 
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(8.7) 


We write the drift function 

by (5.9) with K® given by (7.5), and the transfer function 

T° = - X* <; ez > = - |« - 1 rV) R°(k') , (8.8) 

C ® 

by (5.16b) with X* given by (7.4). Here < t 02 > *= ! ;R° . 

Note 

dK°/dk = 2 B 2 w # -‘ F 22 (k) (8.9a) 

and 

dT°/d(-k) - | “o’ 1 F(k J CR * (k)-R c (k)D 


- k 2 F(k)R* (k) 


(8.9b) 


because r°<< r' in the region of 6mall k where the blocking occurs. 


If we confine the transfer to the spectral component F J2 (k) 
along which the drift operates, the spectral balance (8.6) takes 
the form 


R (k) = <(V u 2 ) > 



( 8 . 10 ) 


after simplification by (8.9a) and (8.9b). The integral equation is 
solved, giving the spectrum 

F 22 (k) = | B 2 k- 5 (8.11) 

and the intensity 

<u 2 2 > = y B 2 k -1 * . (8.12) 
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The eddy coefficient is 


*' »w o “ l e 2 k- 2 . (8.13) 

by (7.4). 

Formula (8.12) indicates that the drift turbulence is limited to a 

length 7 v ~ 1 ( <u ' 2 > /B 2 )^ 

^ 2 


8.3. Transition-range 


A spectral gap lies between the inertial-range and the drift- 
range. It describes the transition from the direct cascade at rate 
the -inverse cascade A and is governed by the balance 

between the two cascades in the form 


This equation governs the enstrophy transfer in a background medium 

having a transport property by (8.13). By writing 

< ; 02 > = JsR°= [ k dkV 2 F(k’) , we solve Eq. (8.14) to obtain 
J 0 

F(k) = 2(u) e e c /B 2 ) k " 1 . (8.15) 

By a comparison between the two spectra (8.11) and (8.15), the 
critical wavenumber that separates the transition-range from the 
drift-range is found to be 

K - (4/3) % 6 (o) e r )-* . (8.16) 

C 0 £ 

The critical waven umb er also determines the passage from isotropy 
to anisotropy. 
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9. Summary and discussions 


The Navier-Stokes equation with the Coriolis force is trans- 
formed into a master equation to describe the microdynamic state of 
turbulence. The group-kinetic method describes the three groups of 
transport processes (spectral evolution, ■ eddy coefficients and 
relaxation) . In the macrogroup dynamics, we derive the kinetic equa- 
tion of turbulence (3.6), the equations of evolution for the vor- 
ticity (5.3) and the enstrophy (5.4), and their governing transport 
functions of transfer T® and drift K° by (5.16) and (5.9). In the 
microgroup dynamics, we derive the eddy coefficients consisting of 

I 

the eddy viscosity K for direct cascade, the eddy coefficient of 
amplification x' for inverse cascade, and the diffusivity K® for 
drift turbulence. These coefficients are calculated from the colli- 
sion integral of the kinetic equation. The subdynamics treats the 
mechanism of relaxation by which the eddy coefficients approach 
their equilibrium. The closure of turbulence is obtained by the 
loss of memory. The subdynamics consists of a Fokker-Planck equa- 
tion (6.1) for path transition and a system of two Langevin equa- 
tions of turbulence (6.5a) and (6.5b). 

All the eddy coefficients K f X and K® are found in (7.3) 
(7.5) to be inversely proportional to rotation u>°. The dependence 
of k' on w {j " 1 has been suggested by Blackadar (1962) empirically. 

The spectral structure of Rossby wave turbulence is investigated 
for the inertial-range and the drift-range as characterized by a 
direct cascade and an inverse cascaide, respectively. The spectral 
gap as formed by the transition between the two cascades is also 
examined. The corresponding spectral laws k* 3 , k“ 5 and k _1 are found 
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in (8.3), (8.11) and (8.15). 


The spectrum k" J in (8.3) and the spectrum k -3 of geostrophic tur- 
bulence (6.8) have separate parameters and should not be confoun- 
ded. This difference cannot be revealed by a theory from the vortex 
model (1.4) . 

It is to be remarked that by the cancellation of and the sub- 
sequent constraint to quasi-geostrophic turbulence in the 
traditional vortex model (1.4), the small eddies cannot organize 
the <d° -dependent cascades to analytically treat the inertial-range 
and the transition range. The drift function cannot find a cascade 
for balance to derive the drift-range. This explains why no theore- 
tical treatment has* succeeded in attacking the Rossby wave tur- 
bulence. 
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GROUP-KINETIC METHOD OF TURBULENCE 
By C. M. TCHEN 

Department of Mechanical Engineering, City College 

New York, NY 10031 

The inhomogeneous Navier-Stokes equation is transformed into a 
master equation in the phase space for the description of the 
microdynamical state of turbulence. The elementary interaction 
among the fluid elements is represented by the pressure-gradient. 
For the three transport processes of spectral evolution, eddy 
viscosity and relaxation, a macrogroup dynamics, a microgroup 
dynamics and a subdynamics are developed. The kinetic equation of 
turbulence is derived. The eddy viscosity is calculated from the 
"eddy collision". The simulation of the subcluster of many high- 
order distributions by an effective fluid medium causes a loss of 
memory. This mechanism of relaxation for the approach of eddy 
viscosity to equilibrium solves the closure problem. The sub- 
dynamics consists of a Fokker-Planck equation of path transition 
and a Langevin equation of turbulence for the determination of the 
orbit function and the path diffusivity. Although the Fourier modes 
interact indefinitely, their statistical quantities interact by 
nearest neighbour-groups. By this property the high-order correla- 
tions enter in the form of a product of pair correlations without 
invoking the hypothesis of normality or relying on perturbation ex- 
pansion and subsequent summation. 

The spectral results are found as follows: (i) In the inertial 

range of the Navier-Stokes turbulence, the energy spectrum is 
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F(k) = 1.650 e 2 / 3 k~ 5 / 3 , and the pressure variance spectrum is 

S(k)=0.907p 2 e 4 / 3 k -7 / 3 , where e is the energy dissipation rate, and p 
is the fluid density, (ii) In the convective-inertial range of heat 

*1/5 *5/3 

diffusion, the temperature variance spectrum is (k)=0. ' k > , 
where e Q is the rate of dissipation of temperature variance. 


1. Introduction 

The derivation of the energy spectrum has attracted the attention 
of many authors. Kolroogoroff (1941) presented a dimensional ar- 
gument, and Heisenberg (1948) proposed an energy transfer from the 
large eddies that present a macro-gradient toward the small eddies 
that organize an eddy viscosity. A particularly appealing method 
has been the renormalization. The group-renormalization transforma- 
tion has succeeded in finding a correction to the 5/3 spectrum from 
intermittencies of turbulence but not in deriving the 5/3 spectrum 
itself (Grossmann § Schnedler, 1977) . This correction is neverthe- 
less too small to be verified by direct measurements. The method of 
renormalization perturbation expansion has enabled Kraichnan (1959, 
1977) to develop his DIA method and generalizations. A systematic 
method of developing the expansion was given in a general form by 
Martin, Siggia 5 Rose (1972), and by Dubois 5 Espedal (1978). 
However, the following difficuties remain: The necessity of the 
molecular viscosity and an artificial driving force as the basis of 
perturbation expansion, the uncertainty of convergence, the 
validity of the lowest-order expansion which amounts to a 
quasilinear theory, and the hypothesis of normality for factorizing 
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a quadruple correlation. 


In order to alleviate these difficulties and clarifying the 
physical roles played by the mathematical operators and functions 
(e.g. evolution operator, collision operator, diffision operator 
and orbit functions) , we develop a kinetic method based upon the 
transformation of the Navier-Stokes equation into the kinetic equa- 
tion of turbulence. We calculate the eddy viscosity from the eddy 
collision-integral, not by perturbation as was done in the Chapman- 
Enskog method for the derivation of the molecular viscosity from 
the Boltzmann equation (Chapman $ Cowling, 1939) . With the know- 
ledge of the eddy viscosity, we derive the spectral structure. 

* 

As the scope of research, we let the pressure-gradient define the 
elementary interaction between fluid elements. The inhomogeneous 
partial differential equation of Navier-Stokes is then transformed 
into a homogeneous master equation in the phase space. The master 
equation has lesser nonlinearity. 

The decomposition into a macro-distribution f°, a micro- 
distribution f’ and a sub-distribution f" represents the three tran- 
sport processes of spectral evolution, eddy viscosity and relaxa- 
tion (Section 2). This decomposition is superior to the reduction 
of the N-body distribution function into singlet-, binary-, and 
triplet-distribution functions in many-body statistical mechanics 
(Bogoliubov, 1962) , because the singlet-distribution f° in the 
group form suffices for describing the spectral function, and the 
groups interact as nearest-neighbours statistically, 

For analysing the three transport processes, we develop a macro- 
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0 1 

group dynamics for f a microgroup dynamics for f and a sub- 
dynamics of relaxation. The kinetic equation of turbulence for * 
is derived in Section 3, and the eddy viscosity is found in Section 
4. The eddy viscosity approaches its equilibrium by relaxation for 
which a subdynamics is developed and consists of a Fokker-Planck 
equation of transition and a Langevin equation of turbulence. The 
former equation governs the path fluctuations and determines the 
orbit function, and the latter finds a path diffusivity (Section 
5) . In sum, the macrogroup dynamics determines the spectral evolu- 
tion from the kinetic equation, the microgroup dynamics and the sub- 
dynamics produce a system of integral equations for the eddy 
viscosity and the path diffusivity. The conversion of pressure 
field fluctuations into velocity fluctuations is made in Section 6. 
In Section 7, we derive the energy spectrum k' 5 / 3 , the temperature 
variance spectrum k' 5 / 3 , and the pressure variance spectrum k -7 / 3 . 
A summary and a discussion are presented in Section 8. 


2. Group-kinetic method 

Fluid turbulence is described by the Navier-Stokes equation 

(S + G-V-vV 2 )G = E, 3,23/31 (2.1) 

/v 

for the fluid velocity u in a driving field i =-^Vp . Here P is 
the pressure, p is the density, and v is the molecular viscosity. 
The incompressible fluid satisfied the equation of continuity 

^ /-v 

V-u=0 . An instantaneously fluctuating function is denoted by ( ) . 

W A A 

For raising the dimensionality we introduce f(t,x,v) * 6Cv-u(t,x)!l and 

v -«». MM. ^ 

transform the inhomogeneous Navier-Stokes equation into the 
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homogeneous master equation 


O + L) f = 0. 


( 2 . 2 ) 


The normalization is 


L = v-V-vV 2 + E • 9 


with 


dv £ = l, and the differential operator is 
9E9/9v The master equation describes 


the microdynamic state of turbulence. A scaling is necessary for 
statistical treatment. We decompose the instantaneously flue- 

A . • ys 

tuating function f = f + f into an ensemble average f =< f > , and a 
fluctuation f . We represent the three transport processes of spec- 
tral evolution, eddy viscosity and relaxation by a macrogroup f 0 , a 
microgroup f' = f* + f" and a subgroup f" = f^+ + respec- 
tively. The operators a, a , A 0 , A*, a" can be used. The groups 

are in the order of decreasing coherence with durations of correla- 
tion t>T° >t' >t" . 

c c c 


Since the spectrum evolves in a medium of eddy viscosity, and the 
latter approaches its equilibrium by relaxation, the three tran- 
sport processes, together with the three groups that represent 
them, interact by nearest-neighbours. This property in association 
with separated coherence constitutes the "nearest neighbour group- 
interaction". Consequently, a quadruple correlation appears in the 
form of a product of two binary correlations as 


< E 1 (t)E 1 (t-T)> <f° (t) f° (t-T) > = <E'(t)E'(t-T)><f 0 (t)f°(t-T)> , (2.3) 

and not in the form < E(t) E(t-x)f (t)f(t-t) > that requires the 

V 

hypothesis of normality for factorization. 

There is a certain analogy between the decomposition of the 
distribution f into groups f°,f, f" and the reduction of the 
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many-point distribution F J2>#iN into the singlet- distribution F l# 
the pair- distribution F 2 , and the triplet-distribution F 3 
(Bogoliubov, 1962) . Both the decomposition and the reduction repre- 
sent the three transport processes of evolution, transport property 
and relaxation that interact as nearest neighbours. However, a dif- 
ference exists in the determination of the spectrum: a kinetic 

equation for f° suffices in the first case, and a system of two 
kinetic equations for F and F 2 is required in the second case. 


The instantaneously fluctuating groups can be transformed into 
Fourier modes with overlapping wavenumbers, but the statistical 
groups are separated by adjacent wavenumber limits 
O^k'^k, k<k"<®, k' r '<k" , <«> as demarcation. Thus 


< E°E°> 


dk' _X<E°(k') E°(-k') > 


< E'E'> {} = jdk" .X<E (k") E (-k") > {} 


<E"E"> {} = 


dk"’ J«E"(k’")E"(-k’")>{} 


l 


It is understood that the volume integrals with respect to 
dk', dk", dk'" extend from - 00 to + 03 and that k separates the 
macro-intensity and the micro-intensity, while other wavenumbers 
are dummy variables of integration. The Fourier decomposition is 

truncated with a truncation factor . 
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3. Kinetic equation of turbulence 

The roaster equation is scaled by A 0 and A 1 to give the transport 
equations 


(3 t +A°L)f° = - L°f +C 1 {f°}, C 1 {f°}=-A°L 1 f 1 (3.1) 

(3 t +A 1 L)f l = - L° (f+f°) + c"( fl } » C"{f*}= - A 1 l"f" . (3.2) 


The coefficients of "eddy collision" C 1 (t,x,v){}, c"(t,x,v){) are 

r - ~~ 

integral operators. By random encounters, the many distributions 

in the cluster f"= f^ 2 -* + f( 3 )+-«- become homogenized and lose their 
individuality in v-dependence . As a result, the cluster can be 

simulated by a frictional medium so that the coefficient of colli- 
sion c"(t, x ,v) { } =C"(t,x) ceases to be an integral operator. 

By (2.3) we can write c Ht.x.v) = c'(t,x,v) and L 1 (f+f°) = L'f°, 
reducing (3.1) and (3.2) in the closed form: 

(3 t +A°L)f° = - L°f + C' { f 0 } , c'{f°} ■ - A°L*f ' (3.3) 

n3 t +A'L(t,x,v) - c'tt.xjH' = -L'f° (3.4) 

*** 9m 

The homogeneous equation 

O t +A , L-C")U(t,t’)= 0, with U(t,t ) = 1 (3.5) 

governs the evolution operator U(t,t ) by means of which we in- 
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tegrate (3.6) to get 


ft 


f = - 


Jo 


dx A U(t,t-x) L (t -x) f° (t-x) 


and find the eddy collision 


(3.6) 


c'{f° (t-x) } = A 0 


ft 


dx L (t , t-x) A U(t,t-x)L (t-x)f°(t-x) 


= 3-D • 3 { f 0 (t-x) } • 


The diffusivity 


(3.7) 


D « 


t = oo 


dx < E ' (t) A'U(t,t-x)E' (t-x) > 


(3.8) 


is adiabatic by the large-time diffusion on account of the separa- 
tion of groups e' and f° . By the knowledge of the collision 
operator (3.7), Eq . (3.3) is called the kinetic equation of tur- 
bulence. 

By taking the moment of (3.3) we obtain the Navier-Stokes equa- 
tion 


(3 t +A°u-V-vV 2 )u° = E°+ J° 


(3.9) 


with the added "eddy damping" 


J° = I dv v C 1 (t ,x,v) (f° (t -x) } . 


(3.10) 


By multiphying by u° and averaging, we find the equation for the 
spectral evolution in quasi-homogeneous turbulence 


1 9 < u 02 > 

2 t 


- v<(Vj u?) 2 > - T° 


(3.11) 
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with the transfer function 



(3.12) 


4. Kinetic origin of the eddy viscosity 


4-1. Transformation of the evolution operator into orbit functions 


At the instant of time t-x the Lagrangian field 
A'u(t,t-T)E' (t-x)= a'e' Ct-T,x(t-x)D takes the position 

x(t-T) = x - vx - £(x) along the perturbed trajectory that 

passes through the phase point x, v at the time t . Here £(t) is 
the path length traveled in the interval of time x . Thus the 
Lagrangian correlations are 


< U°(t)f Q (t -T )> 


dk' X < u°(k , )f°(-k , .v)> h v (x,k', v)h^(x ,k ') 


(4.1) 


< e' ( t)E'(t-x)> = j dk" X < E , (k")E’(-k")> h^x.k", vjh^x.k") 


f.. - 




(4.2) 


= dk’”, 


< E"(t)E"(t-x)> = |dk'”X < E"(k"')E" , (-k'") > h^T.k") 


(4.3) 


The field correlations can be integrated to form the adiabatic dif- 
fusivities: 


d' ( t,x,v) 



dk" X< E’(k") E'(-k")> (x ,k", v) h^(x ,k") 


(4.4) 
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D (t,x) = 


dx j dk" *< E"(k") E"(-k")> hJ(T,k") 


(4.5) 


The orbit functions are: 


h (T,k’,v) = exp(-ik"*v T ) 
V tm, 


(4.6) 


hjl (T,k'') = EexpC”(k”)T 3 <exp-ik”-r(T)>, hJ(T.k> <exp - ik ,, -£ ,, (T)> < 4 ‘ 7a ' b) 


h (T,k") = <exp- ik"* (fi,° + £. 1 )> = exp(- j k" u 2 t 2 ) 
m - - ^ 


(4.8) 


Here h y governs the streaming, and govern the relaxation by 

£" and h represents a meandering by 8° and 8 with an effective 
^ m ,a ~ 

velocity u°. The meandering without relaxation appears in the lowest- 
order perturbation expansion in quasilinear theory. The strong tur- 
bulence is predominantly governed by streaming and relaxation. 


4.2. Transfer function 

By substituting (4.1) and (4.4) in (3.12), we get the transfer 
function 


fj>0 _ ^ 


dv v. 
i 


dk"dk’ ^<E , (k")E'(-k")> h^Cr.k'^h^Ci.k') 


X | dx M° (T,k",k' ,v) { X <u? (k')f(-k ,v)>} 


(4.9) 


with the memory operator 


M°{ 


a h v (T,k", v) a^CT.k'.v) 


(4.10) 
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In differentiations 3. 9. we find terms that cancel upon integra- 

3 3 

tions with respect to dv dk" dk'. The remaining terms reduce the 
memory operator to: 


(i) transport in a macro-gradient 


M°{ } = h (T,k",v) 9 2 h ( x,k* ,v){}, i-im M°{} = - k ' 2 t 2 
V - - V ~ ~ v=0 

(ii) transport without gradient 


(4.11) 


M° { } = - 9 2 h ( T , k", v) h (T,k'v){), Zim M°{) = k" 2 T 2 . (4.12) 

v v - * v=0 

r 

Note that (4.11) and (4.12) govern the memory transmission by 
v-dependence through h . The loss of memory at the limit v = 0 can- 
cels the role of operator and reduces the transfer function into 


k' r° 


for direct cascade 


o _ 


-A.’ < u 02 > for inverse cascade 


(4.13a) 


(4.13b) 


by (4.11) and (4.12), respectively. In the direct cascade, energy 
is transferred from the large eddies that present a macro-velocity 
gradient such that 


R° 


<(V.u?) 2 > 
3 i 



dk k 


F (k ') 


towards the small eddies that organize an eddy viscosity 


(4.14) 
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(4.15) 


f r°° 

K' - tr k' - j dk’tr X<E’(k') E'(-k') > J^dx T 2 h £ (T,k') 


Here F(k) is the energy spectrum. In the inverse cascade, the 
small eddies organize an eddy amplification at a rate 


A* = J dk" k" 2 tr a < e' (k")E_' (-k") > dx T 2 h £ (x,k") (4.16) 

and feed energy into the large eddies. 

The approximation h £ (x,k")h £ (x,k') = h £ (T,k”) by k'« k" has been 
made. The moment [ dv v f° = u° has been taken, and other details 
of calculation have been omitted. In the following we confine our- 
selves to the transport by K ' in a direct cascade only. 

For the sake of abbreviation, we write 


dx t 2 h p ( t , k ' ) 
Jo * ~ 



= 2y' " 3 , Y’(k’) = c ’ (k 1 ) ♦ Y £ (k') 

r°° 

dx<exp-ik"A"(T)> y £ _1 , Y £ (k") = k" 2 K £ . 


(4 . 17a, b) 


(4 . 18a , b) 


The subdynamics for the analytical derivation of (4.18b) is given 
in Section 5. 


The relaxation times y' 1 and y £ 1 
equilibrium for the eddy viscosity 


govern the 


approach to 


k’ = 2< e'e’> {Y' ~ 3 } 
and the path diffusivity 


(4.19) 
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(4.20) 


<A"(t) £"(t)> = < e"e" > {y'^ 3 } , 

Y = 00 ~ "* 


respectively. 


4 . 3 Eddy coefficient 

Recall that upon simulating f" by an effective fluid medium, we 
have reduced the collision integral C"(t,x,v){} = C"(t,x) into a col- 
lision coefficient C"(t,x). If the same reduction is made for 
c'(t,x,v) =C'(t,x) , we reduce the eddy damping (3.10) into 

J° = c'(t,x) u°(t,x) and the transfer function (3.12) into 


T° = - < u°* c'(t,x) u° > • 

By comparing with (4.13a) we identify 
C'(t,x) = k'v 2 
and similarly 

C"(t,x) = K"V 2 , or C"(k) = -k 2 K" . 

Hence we find 


(4.21) 


(4.22) 


(4.23) 


Y ’(k’) = k ' 2 (k' 
from (4.17b), and 


+ " ( p t + ^ 

transform (4.19) 


into 


(4.24) 


K' = 2(P t +l)’ 3 tr <E 'e'> P t = K’/K^ . (4.25a,b) 

The path diffusivity governs the orbit functions in (4.17a) and 
(4.18a), and controls the approach of the eddy viscosity (4.25a) to 


13 



equilibrium. By a subdynamics in Section 5, we shall analyse this 
transport coefficient. 

5. Subdynamics of relaxation 

The orbit function < exp -ik"* £"(x) > and the path diffusivity 
govern the relaxation. For their determination we develop a 
subdynamics consisting of a Fokker-Planck equation of transition 
and a Langevin equation of turbulence. These equations are in the 
configuration space by definition of relaxation. 

5.1. Fokker-Planck equation of transition 

The probability that a path-length is made between £ and £ + d£ in 
a time interval t along a trajectory that passes through the point 
x at the time instant t is 


p(t,x;t-T, x-£)d£ or briefly 


p(T,£)d£ 


in quasi-stationary processes. The 
trograde transition is governed by 
(Tchen, 1944) 


probability density of re- 
the Fokker-Planck equation 


frpt T -i> * K * (5-D 

** » 

and satisfies the condition of normalization 

jd£ p(x,£) = 1 . ( 5 - 2 > 

The streaming is absent in relaxation. The path diffusivity as 

defined by (4.20) is adiabatic, i.e. independent of x and £ . 
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By Fourier decomposition, we transform Eq. (5.1) into 


$fP (T.V") = - k" k": Kj p(T,k") 

with the solution 


(5.3) 


p(x,k") = — - — exp[-k" k": K. t] . 
(2tt) 3 ~ ~ = 


(5.4) 


The coefficient of integration is determined by the condition 
(5.2). Hence we obtain 


f -ik •£ 


< exp -ik"*£ (i) > = 


lx e 


P(T.A) 


= exp(-Y^ T) , 


" n, " 

and = k . 


(5 . 5a,b) 


The path diffusivity and the relaxation time as determined by 
(5.5a,b) have been introduced in connection with (4.17) and (4.18). 


5.2. Langevin equation of turbulence 
The Langevin equation of turbulence 


du"(k") 


dt 


" 'Vi r" /-l '\ 

+ u (k ) = E (k ) 


(5.6) 


governs the relaxation in k"- space. It has a certain analogy with 
the Langevin equation for Brownian motion. 


By two successive integration we get 
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1! 

Z = 


dt 


^ dt"expC-Y"(t '-t") D e" ( t") 
0 x ~ 


ft If ft I If I II II II 

= j o dt j M dt exp[-Y £ (t -t )] E (t ) 


ft 


■ Y c 


dt" E"(t")(l - expC-y’JCt-t")]) 
o — *• 


(5.7) 


An interchange of order of integrations has been made. In cal- 
culating k'^ by definition (4.20), we find four double integrals, 
one of which is 

>. i i 


t* — 1 

= Y P 

r t ft 

dt’ 1 dt"< 

|E"(t ')£>") |> 

X 

0 J 0 



.ft . ft' 

. M 11 , 

= 2Y *■ dt dt 

0 0 

< | E (0) E (T)|> 


r 

ft 


= 2y"” 2 

j Q dT(t-T) < | 

E ( 0) E (T)|> . 


(5.8) 


An interchange of order of integrations leads to 


j Zim-£<\z'\t)Z'\t)\> = 


2 


f t-^°° 
0 


dx< |e"(0)E"(t) I > 


(5.9) 


or 


K ; . .>x-> , (5.io) 

by (4.20). Note that the remaining three double integrals decay 
exponentially and do not contribute. Also note that (5.9) is 
written in k"-space, by (5.6), i.e. 


<|A"(t)a"(t)|>«< z'\t ,-k")>, <|E"(0)E n (T)|> = <E’ , (t,k”)E n (t-T,-k’ , )> 
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It is to be remarked that a relationship analogous to (5.10) ex- 
ists in Brownian motion and is called Einstein's formula of fluctua- 
tion dissipation. 

By substituting for d" from (4.5), we obtain 

K Jt “ tr K A = tr < f ^4 * (5.11) 

and confirm (4.20b). 


A comparison between (5.11) and (4.25a) yields the quartic equa- 
tion 

r 

P t (P t + l) 3 = 2 • (5.12) 

The solution determines the turbulent Prandtl number P t =0.545. 

6. Eddy coefficients as functions of energy spectrum 


6 . 1. Equation of state of turbulence 

By the condition of incompressibility V • u = 0 the Navier-Stokes 
equation yields the equation 

V*E = r, or E (k ) = l k k r(k) (5.1) 

with r" = a" VV: u u . It applies to relaxation in equilibrium. 
Upon neglecting the non-stationary effects by meanderning, we decom- 
pose Vu(jc) = Vu°+Vu"(x) into a quasi-homogeneous background V u° and 
a variable sub-gradient Vu"(x), and write 



( 6 . 2 ) 


r" = A (V.u.+V.u.) (V.u: +V.u. ) . 

13 1 3 3 1 31 


By selecting the subgroup, and assuming isotropy, we get 


< r" 2 > = §- R° r", with R° = <(V.u°) 2 >, R" = <(V.u") 2 > 

9 j l' 3 1 

and derive 

< E" 2 > = | R°< u " 2 >, 
by ( 6 . 1 ). 


(6.3) 


(6.4) 


6.2. Eddy coefficients 

r 

By writing (6.4) in the operator form < e"e">{} = |- R°<u" u" > {}, we 
derive 


K £ = | R ° tr < u" u">{y”' 3 } 


(6.5) 


by (5.11), or 


K^k') = |-R°(k) K^k^UY^k")]- 2 } 


( 6 . 6 ) 


by = tr <u"u">{y£ *} from definition (4.20). 

The integral equation ( 6 . 6 ) can be solved by a differentiation 
with respect to k', yielding 


Vk) = (2/9) h R oh k ' 2 


(6.7) 
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It follows 


from (4.25b) and (6.7). The numerical coefficient is 

c„ = (2/9)** P = 0.471, for P = 0.545, 

by (5.12). 

7. Spectral structure 


7.1. Energy spectrum 

r 

Formula (6.8) for the eddy viscosity K helps to determine the 
energy transfer across the spectrum. In the inertial-range, the 
transfer at a constant rate of dissipation e is 

k' R° = e, or (2/9)^ P t R 03/2 k" 2 = e, (7.1) 

from (3.8) and (4.13a). The vorticity function R° is defined by 
(4.14). By isolating R° and differentiating, we obtain the Kol- 
mogoroff law 

F(k) = A e 2 / 3 k -5 / 3 (7<2) 

2 -2/3 

The numerical coefficient is evaluated to be A = J C k 1-^0 , 

by (6.9) . 

7 . 2. Pressure fluctuations 


The pressure fluctuation p _1 p(k) = -ik 2 k-E(k) 


can be cal- 



culated from (6.1), giving 


X < | p (k) | 2 > = p 2 k 2 OC < |E (k) | 2 > 


= fp 2 R°(k) k -2 *<|u"(k)| 2 >, 


(7.3) 


by (6.4). 

By introducing the spectrum S(k) of pressure fluctuations such 


that 


2 j dk" S(k") = jdk" X< |p"(k") | 2 > 


we can write (7.3) in 


the spectral form 


S(k) = |p 2 R° k" 2 F(k) 


= C e 


2 r 4 / 3 k ' 7 / 3 


(7.4) 


with C = A 2 = 0.907 . The use of (7.2) has been made. 


7.3. Thermal spectrum 


The evolution of a passive scalar, e. q. the temperature fluctuation 


'V 


6, is governed by the homogeneous diffusion equation 


O t + 'k u • V - KV 2 )fr = 0 


(7.5) 


in isotropic turbulence. The heat transfer in the convective- 
inertial range is described by 


kI R° a = e Q or (2/9)^ R° Q = e Q , 


v l 0 0 


(7.6) 


20 


where e fl is the dissipation rate by the molecular diffusivity <, 


and R° = <(V9 0 ] 2 >= 2 


rk 


dk' k' 2 F(k') is the variance of the temperature 


gradient. 


By the same procedure as was used for deriving (7.2) from (7.1), 
we obtain the solution of (7.6) as 


F 0 (k) = B e 0 e' 1 / 3 k' s/3 . (7.7) 

The numerical coefficient is found to be B = 0.899 in consistency 

with the condition B = P A. 


8. Summary and discussion 

By letting the pressure gradient represent the elementary interac- 
tion among fluid elements and by raising the Navier-Stokes equation 
to higher dimensionality, we obtain the master equation for the de- 
scription of the microdynamical state of turbulence. Being a 
homogeneous equation, the master equation is suitable for the ex- 
plicit formulation of the eddy transport in the gradient form 
(direct cascade) and in the form of an eddy damping (inverse 
cascade) . By the macrogroup dynamics and the microgroup dynamics, 
the kinetic equation is derived, and the eddy viscosity is found to 
depend on the orbit functions of streaming and relaxation, in the 
roles of memory-transmission and memory-loss, respectively. The 
loss of memory leads to closure and irreversibility. 

The subdynamics consists of a Fokker— Planck equation of transi- 
tion and a Langevin equation of turbulence, and governs the path 
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perturbation along an accelerated trajectory and the path dif- 
fusivity . 

The eddy viscosity and the path diffusivity are found to be 
governed by a system of integral equations. The solutions deter- 
mine the energy spectrum F(k) = 1.650 e 2 / 3 k 5/ ^ 3 , the pressure 
variance spectrum S (k)=0.907p 2 e 4 / 3 k -7 / 3 in the inertial range, and the 
temperature variance spectrum F Q (k) = 0.899 e Q e ly ^ 3 k 5 / 3 in the 
convection-inertial range. 

If the orbit function is degenerated into a non-adiabatic diffu- 
sion as the lowest order perturbation expansion in the quasil inear 
theory, a meandering at an effective velocity u q is found. The 

* u , -3/2 

parameter eu ( yields a spectrum F(k) % (eu ( ) k reproducing the 

result of DIA. 

Since the groups represent the transport processes of spectral 
evolution, eddy viscosity and relaxation, and are analysed by a 
macrogroup dynamics, a microgroup dynamics and a subdynamics se- 
parately, the high-order correlations appear in the form of group- 
interactions <u' u' > <f° f°> by our group-kinetic theory and not in 
the form of a quadruple correlation <uuf f > as would be by a per- 
turbation theory. The quadruple correlation requires the 
hypothesis of normality for factorization and closure. 

9 Acknowledgment 

This work was sponsored by the National Aeronautics and Space Ad- 
ministration under contract NAS 8-35153. 


22 



10. REFERENCES 


BOGOLIUBOV, N. N. 1962 Problems of a dynamical theory in statistical 
physics, in: Studies in Statistical Mechanics . I ed. J de Boer 
§ G. E. Uhlenbeck, North-Holland, Amsterdam. 

CHAPMAN, S. 5 COWLING, T. G. 1939 The Mathematical Theory of Non- 
Uniform Gases , Cambridge University Press. 

DUBOIS, D. F. 5 ESPEDAL, M. 1978 Direct interaction approximation 
and plasma turbulent theory. Plasma Phys . 20 , 1209-1239. 

GROSSMAN, S. $ SCHNEDLER, E. 1977 Fluctuation correlation of the 
turbulent spectrum. Z. Physik B 26 , 307-317. 

HEISENBERG, W. 1948 Zur Statistichen Theorie der Turbulenz, Zeit . 
Physik 124 . 628-657. 

KOLMOGOROFF, A. N. 1941 The local structure of turbulence in 
incompressible viscous fluid for very large Reynolds number. 

C. R. Acad. Sci . URSS , 30, 299-303. 

KRAICHNAN, R. H. 1959 The structure of isotropic turbulence at very 
high Reynolds numbers. J. Fluid Mech. , 5, 497-513. 

KRAICHNAN, R. H. 1977 Eulerian and Lagrangian renormalization in 
turbulent theory. J. Fluid Mech. , 83, 349-374. 

MARTIN, P. C., SIGGIA, E. D. 5 ROSE, H. A. 1973 Statistical dynamics 
of classical systems. Phys. Rev. , A8 , 423-437. 

TCHEN, C. M. 1944 Enige wiskundige betrekkingen welke een rol spelen 
in dif fusie-problemen (Some meathematical relations which play a 
role in diffusion problems) . Versl. Ned. Akad. Wetensch. Afd. 
Natuurkunde, 53, 400—410. 


23 



